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We propose a possible charge fluctuation effect expected in layered superconducting ma- 
^0 , terials. In the multireference density functional theory, relevant fluctuation channels for the 

. Josephson coupling between superconducting layers include the interlayer pair hopping de- 

rived from the Coulomb repulsion. When interlayer single-electron tunneling processes are 
j-H , irrelevant in the Kohn-Sham electronic band structure calculation, the two-body effective 

O , interactions stabilize a superconducting phase. This state is also regarded as a valence-bond 

^ , solid in a bulk electronic state. The hidden order parameters coexist with the superconduct- 

Oh, ing order parameter when the charging effect of a layer is comparable to the pair hopping. 

c/3 , Relevant material structures favorable for the pair-hopping mechanism are discussed. 

^ ' KEYWORDS: superconductivity, iron pnictide, cuprate, density functional theory, Haldane gap 

^ ■ 
^ : 

O . !• Introduction 

O . 

The discovery of iron-pnictide superconductors gave us an interesting ground for testing 
' theoretical approaches to analyze superconductivity. The first records of jump in the supercon- 

i ducting transition temperature in iron pnictides was observed in fluorine-doped LaFeAsO^'^^ 

. when the discovery of fluorine-doped LaFePO^-* fuelled a search for chemical trends in series 

■ of superconductors. Similar tests are also necessary to develop the theory of high-temperature 
O ' superconductivity. 

I In the present standard of theoretical approaches to analyze superconductivity in materials 

. science, we may apply the strong-coupling theory of superconductivity starting from estima- 

^ ■ tion of electron-phonon coupling constants e.g. by using the density- functional perturbation 

' theory (DFPT).^-* The applications of DFPT in high-temperature element superconductors^' 

suggest that the technique enables us to obtain a reasonable estimation of the transition tem- 
perature Tc when the essential pairing mechanism is the electron-phonon-interaction-mediated 
stabilization of superconductivity. The method may also predict how to enhance the Tc of com- 
pound superconductors. As an example of its application, I and my coworkers have shown that 
a possible positive jump in in CaSi2 is expected at a high pressure when a structural phase 
transition to the AIB2 structure takes place. ''^ 

Recently, Jishi and Alahyaei have used this method to analyze the iron-based superconduc- 
tors LiFeAs and NaFeAs.^-* These 111 compounds may be ideal arsenides for the application 
of the strong-coupling theory, since clear superconducting transition temperatures of 18 K for 
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LiFeAs^'^°) and 9K for NaFcAs^^) are reported under ambient condition without doping. An- 
other experimental report on NaFeAs suggests a higher transition temperature of above 12 K 
without a clear indication of any coexisting magnetic ordering. -"^^^ Jishi and Alyahyaei unfor- 
tunately failed to obtain the transition temperature observed in experiments and concluded 
that the theory indicates a Tc far below 1 K owing to insufficient electron-phonon coupling 
constants. 

If one looks at a list of related materials/'^) he or she might be motivated to perform 
another theoretical test in order to understand a chemical trend. The technique was pplied 
to NaFeAs by A. Nakanishi and he considered a chemical trend by studying other two sup- 
posed material structures: NaCoAs and NaNiAs. Using the calculation results, Kusakabe and 
Nakanishi studied the electronic structures of these compounds to derive a chemical trend. ""^^^ 
A key to understand the trend could be found in the direction to substitute Fe with other 
transition-metal elements, since electronic bands around the Fermi level are known to origi- 
nate from iron 3d orbitals in Fe2As2 layers. However, note that there are only a few examples 
of nickel arsenides in the 1111 structure showing a finite of approximately 3K}^'^^^ For- 
tunately, Nakanishi succeeded in observing the theoretical stability of the hypothetical com- 
pounds NaCoAs and NaNiAs, as well as NaFeAs, using first-principles structural optimization 
techniques^'^^ with the generalized gradient approximation (GGA).^^^ He applied the strong- 
coupling theory and estimated the transition temperature within a standard approach. By 
means of first-principles lattice dynamics, the superconducting transition temperature was 
estimated to be approximately 0.034 K for NaFeAs, assuming that the system remained non- 
magnetic. An interesting finding was the increase in Tc in the order of NaFeAs (0.034 K), 
NaCoAs (0.127 K), and NaNiAs (3.573 K), although there might remain systematic errors due 
to the limitation in the computation. However, this possible chemical trend is in inverse re- 
lation to the experimental finding of Tc =9K or Tc >12K for NaFeAs and the Tc for Co and 
Ni compounds not higher than that for this iron arsenide. 

As for the Kohn-Sham band structure on NaFeAs, ^'^^ the essential features around the 
Fermi level are rather similar to those of LiFeAs.^^'^'^^ The effective bands in a nonmagnetic 
solution are almost dispersionless along the c direction at approximately the Fermi level. Since 
electron and hole branches coexist and since DOS is reduced at approximately the Fermi level, 
we can conclude that the band structure has a semimetallic nature. That is, it shows strong 
two-dimensionality. Differently from LiFeAs, NaFeAs has two hole pockets. LiFeAs has three 
hole pockets on the two-dimensional plane including F-X-M lines in the first Brillouin zone. 
In addition, the two-dimensionality of the Fermi surface is much better for NaFeAs, and all 
of the pockets of NaFeAs except for one hole pocket around the T point are rodlike. 

A crude realization of the rigid band picture among three band structures for NaMAs 
(M=Fe, Co, Ni) has been reported. The correspondence of the major branches is clearly seen. 
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Thus, the band structures of NaCoAs and NaNiAs arc approximately a result of heavy-electron 
doping in the band structure of NaFeAs. In many iron arsenides, we find a semimctallic two- 
dimensional band structure. ^^"^^^ Following this interpretation, we can understand the reason 
why the electronic density of states (DOS) at the Fermi energy increases when one considers 
the DOS's of NaCoAs and NaNiAs in comparison with that of NaFeAs. The Nakanishi data 
suggesting the increase in Tc for NaCoAs and NaNiAs comes from the larger DOS together 
with the presence of more three-dimensional Fermi surfaces for these two compounds than for 
NaFeAs and the assumption of the electron-phonon-interaction-mcdiated superconductivity. 

To go one step further, we can start analyzing the characteristic topology of Fermi sur- 
faces. The chemical trend of the Kohn-Sham band structure in the series of NaFeAs, NaCoAs, 
and NaNiAs in PA/nmn suggests that rodlike two-dimensional Fermi pockets appear only in 
NaFeAs, while more three-dimensional characters for Co and Ni compounds indeed enhance 
Tc if the strong-coupling theory with the electron-phonon coupling is assumed to be applica- 
ble. This result as well as the theoretical data on the above known first-principles electronic 
structure calculations of 111 compounds leads us to a conclusion along the following line. An 
important ingredient tractable in the density functional theory (DFT) is charge fluctuation 
modes. If one of the modes becomes relevant on a two-dimensional Fermi surface, and if the 
fluctuation effect enhances the stability of a superconducting state, the theoretical approach 
starting from the standard Kohn-Sham scheme^^^ would be feasible. 

In this study, we investigate charge fluctuation effects tractable in the multireference 
density functional theory (MR-DFT).^'^^'^'^) The two-dimensional electronic structures of iron 
arsenides found in Kohn-Sham band structure calculations actually suggest a Cooper-pair 
hopping mechanism in layered materials. We derive an effective theory of superconductivity. 
Important point is weak interlayer single-electron hopping processes in iron arsenides. In 
our MR-DFT formalism, charge fluctuation modes are introduced by flxing the Kohn-Sham 
single-particle description as a mean-field limit of the theory. Here, a simple effective Bosonic 
Hamiltonian is derived for layered superconductors. The model tells us that the formation 
of the valence-bond-solid state in a bulk superconducting state enhances the stability of the 
ordered state via the appearance of hidden order parameters. A hypothesis on the stabilized 
superconducting state will be addressed, where minimized charge fluctuation leading to the 
S = 1 effective one-dimensional Heisenberg spin Hamiltonian is required for the most stable 
superconducting state in layered materials compared with other Heisenberg models with S > 
1. This picture is confirmed if we assume that quasi-particles in the layered superconductor 
highly correlate. Discretized quasi-particle spectrum expected in correlated electron systems 
ensures the effective spin Hamiltonian. 

To start discussion about the first-principles simulation method concrete, we address a 
new theory for the correlated electron systems called the density functional variational theory 
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(DFVT). This is a simple variational method that always refers to a self-consistent solution 
given by MR-DFT. Finally, a means of applying the analysis techniques of DFVT to iron 
arsenide, as well as to other layered high-temperature superconductors including cuprates 
and MgB2, will be addressed. 

2. The Multireference Density Functional Theory 

In the standard DFT, the explicit form of the charge fluctuation is given by the energy 
density functional,^^) 

EfiucH = |^'dA|dVdV^^(^'^|(n(r)-n(r))(n(r')-n(r'))|*^) 

Here, n(r) is the electron density, '^^(r) is the electron field operator for the spin a, n(r) = 
V'CT(r)V'(T(r), and is a minimizing wavefunction of the following reduced energy density 
functional 

Fx[n] = min (V|f + AKel*')- (2) 
In the above definition, the kinetic energy operator T is 

where m is the electron mass, and the Coulomb interaction is given in operator form as 

t>ee = ^y"d3rdV^^^5;V'i(r)Vi(r')V'a'(r')V'.(r). (4) 

a,a' 

Although eq. (1) is an exact expression, it is not easy to obtain insight into relevant charge 
fluctuation modes in a solid only by considering the use of this functional form. This is partly 
because the form is written in double integrals with respect to electron positions r and r'. 
However, we can see that the charge fluctuation may occur everywhere in the electron system. 
Once a crucial scattering due to the Coulomb fluctuation occurs between Fe2As2 layers, a pair 
scattering from one layer to another becomes a forward scattering. 

In our multireference density functional theory, we can introduce part of the fluctua- 
tion explicitly by introducing a fluctuation term, (^'|VxJ*), with the self- interaction correc- 
tion in the form 

n ^ 

Since we may have a series of models, we use the notation (*|VxJ^') with index i specifying 
the number of the model. The notation : O : with an operator O denotes the normal ordering 
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with respect to the creation and annihilation operators. The operators i^^"-* and Z^^^^ may be 
given (i) by an expansion formula of the Coulomb operator around a fixed center and (ii) by a 
creation method for the Dirac character for the crystal. Definitions of each operators are given 
in the literature. Following the work,^^) we introduce the notation <l>pLM(r) = gpL{'>^)YLM{G, <p) 
and ^pLM(r) = 9pL{r)YLM{(^-,4') using the spherical harmonics Ilm and a complete set of 
9pl{t) expanding the radial waves. Another function cjpLir) is 

~9pL{r) = f dr'{r')^gpL{r'). (6) 

^ Jo 

In our Coulomb operator expansion formula, the Y and Z operators are given by, 

YpLM = -^J d^rY^ Vi(r) ($pLM(r) + $pLM(r)) Va(r) , (7) 

ZpLM = -^J {^pLM{r) - %LM{r)) Va(r) • (8) 

The nth operators Y-^"^ and zj^"'^ in the ith. model may be given by identifying a parameter 
set p, L, and M as n = {p, L, M). We have other possible expansions using a screened form of 
the Coulomb operator if we utilize DFVT given in the next section. The definitions of eqs. (7) 
and (8) are independent of the Kohn-Sham orbitals, which are used to expand the wave- 
functions and the field operators in the creation and annihilation operators. This point gives 
an advantage to our formalism, because the scattering channels are defined before obtaining 
expressions of the Kohn-Sham orbitals. 

We also use the notation to represent a density associated with a state | ^) as n^i (r) = 
(^'|n(r)|^'). The energy functional of the new extended Kohn-Sham scheme is 

GxJ*] = (*|r + T6f,|*)- min {^'\f + Vx^\^') 

(r) 

+ F[n^] + j (i^ruext(r)n^(r) 

+Er^c[n^] + j (i^ri;ext(r)n^(r) . (9) 
Here, we refer to the universal energy functional F[n\ given by 

F[n] = min(^''|f + 1/ee|^')- (10) 

The definition of -E'rxc['^] is given by eq. (9) itself. The new extended Kohn-Sham model is 
actually an effective many-body system. 

When we let some of X^^^ be finite, if we replace -Erxc[?^] with the GGA energy func- 
tional EGGA[n^\-, the model becomes a correlated Fermion model, which is defined by the 
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approximated energy functional 



r — r 



+£;GGA[n*] + y d^TOext(r)n*(r). (11) 

The two-body scattering process happening in the z-th model is derived from the fluctuation 
term, which is divided into the effective two-body interaction and a counter term as 

— /Tj -t- n.j (.Qunter ■ \^^J 



n ^ ' 



(13) 



n 

-((if^)V^f^-(^fy-(^i"^)|- (14) 
Since we have the exchange-correlation potential for GGA as 

a secular equation is derived by imposing the normalization condition of 1^*) using the La- 
grange multiplier E as, 

= nfm = Em, (16) 

= f + J v^s{r)n{r)d^r . (18) 
The effective single particle potential Veff(r) is given by 

= / y^l'^^''' + '"GGA(r) + i;ext(r) • (19) 
We have a potential problem given as, 

|-^^r + Weff(r)| Xm,k,fe.(r) = £m,k,feXm,k,fe (r), (20) 

in which Bloch orbitals Xm,k,fe^(r) are determined to be normalized and orthonormal in a 
crystal phase. Here, the Kohn-Sham orbital is specified by a two-dimensional wave vector, k, 
another wave vector, kz, along the c axis, and the band index m. 

We do not explicitly write the i dependence on H^. However, it is implicitly depen- 
dent on the fluctuation term (^'|Vx.|*) through the self-consistency on the charge density 
n(r) = n*(r). We may have another definition of a single-particle problem by including the 
nonlocal potential part and/or mean-field part coming from {Y^i^^) and {zj^^^}, i.e. the 1-body 
counterterm "7^ J counter- However, this shift only changes the definition of single-particle or- 
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bitals used to determine the many-body problem. The charge density n,j,(r) determines the 
expectation values {Y-^"^) and (Z^^). Thus, if n,j,(r) is almost unchanged in a self-consistent 
loop, these values also remain unchanged. Even if there are many-body correlation effects 
creating a few meV of gap in the model, the essential features of the GGA band structure 
are unaffected by only the part WJ^,Q^^^gJ.. An essential change may occur via the correlation 

effects appearing as a slight change in n^{y) and a large change in variational energy via "H?. 

ill) 

When X\ = for all n values, we obtain a secular equation of the Kohn-Sham equation 
in GGA given by eq. (20). The ground state of this model is slightly shifted from the final state, 
when we consider a correlated electron system. The obtained single-particle description defines 
the Kohn-Sham band structure. Introducing a proper Fourier transformation for a selected 
set of bands, we immediately obtain the Wannier representation, allowing us to rewrite 
in a second-quantized form of a tight-binding model. 

In the GGA calculation, we may consider the spin density by introducing a spin-dependent 
GGA functional. If we adopt the spin-GGA scheme, the Kohn-Sham equation becomes 
spin-dependent for a magnetic solution. For the formulation of MR-DFT, however, spin- 
independent Kohn-Sham orbitals are very useful for the discussion of magnetism and su- 
perconductivity. This is because the correlation effects are described by the multireference 
variational states, i.e. the multi-Slater determinants, in MR-DFT. Thus our starting point is 
a paramagnetic state obtained by a nonmagnetic GGA calculation. 

3. Pair-Hopping Mechanism 

In the first-principles study of compound superconductors, we may apply the Kohn-Sham 
scheme^^^ of the density functional theory (DFT)^^) as a starting point. Thus, we investigate 
Kohn-Sham band structures, which are known in the literature or found in actual calculation 
results. The important points considered here are the structures of the Fermi surface and the 
dimensionality of the low-energy branches of the band structure. 

The pictures drawn from the study summarized in the last section are as follows. We have 
several examples of iron arsenides, which have two-dimensional Fermi surfaces, within the 
Kohn-Sham scheme. Details of the Fermi surfaces actually depend on the type of material. 
However, in general, the GGA band structure suggests a picture of a stack of Fe2As2 layers 
weakly coupled by single-particle tunneling at the Fermi level. 

If we expect Fermi instability due to the remaining Coulomb fluctuation modes lost in 
a mean-field approach, or if relevant fluctuation only shifts the electronic structure around 
the Fermi level by opening a gap of a few tens meV, only a slight change in the total charge 
density due to appearance of secondary order parameters is expected. Then, a picture of the 
layered two-dimensional Fermi gas system should remain as a starting one-body mean-fleld 
state even in the final many-body solution aflFected by correlation eflFects. 

In this picture, we need to incorporate effective two-body interactions both within a layer 
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and between layers. Here, note that there is no interband single-particle hopping process, since 
the single-particle part is diagonal in the band index. A step to maintain self-consistency is 
necessary after solving the many-body problem given by the effective two-body repulsive 
interactions, since the charge redistribution might or might not modify the effective single- 
particle excitation spectrum at approximately the Fermi level. 

The stacking of two-dimensional systems in one direction forms a layered bulk system. The 
system becomes a correlated electron system, because the single-electron tunneling process is 
reduced owing to the two-dimensionality and because possible charge fluctuations in a layer 
should be suppressed also owing to the localized nature of Fe 3d orbitals in the superconductor. 
This picture is tractable by constructing a one-dimensional Wannier representation from the 
Kohn-Sham orbitals. The band dispersions around the Fermi level are almost flat along the 
c-direction, so that the one-dimensional Wannier representation is natural. 

The motion of electrons in a layer is described by a two-dimensional electron gas model. 
Here, a two-dimensional wave vector k and a band index m, or a combined index j = (m, k), 
are used to specify each Wannier state in the Z-th layer. In the two-dimensional system with 
multicolored Fermion quasi-particles in the layer, a reduced model might be well-described by 
a known two-dimensional model. In such a model, the effective electron-electron interac- 
tion parameters would specify intralayer scattering processes. A superconducting fluctuation 
may occur, if we consider the two-body scattering processes due to charge or spin fluctua- 
tions, which is not explicitly counted in the GGA calculation. The electron-phonon interaction 
may contribute to the stability of the superconducting fluctuation in the layer. For our next 
discussion, however, intralayer effective attractions can originate from any mechanism, as far 
as it effectively supports the formation of the precursors of the bulk superconducting order 
parameter. 

Here, the important driving force derived in this study is the interlayer pair-hopping pro- 
cesses between the layers. When the Kohn-Sham orbitals are determined using the effective 
exchange-correlation potential, Va:c(r), these charge fluctuation modes are not included explic- 
itly. In MR-DFT, the modes are explicitly introduced using two-body operators. The two-body 
processes are given by the repulsive nature of the Coulomb electron-electron scattering. If we 
have any effective attractive interaction at a position between the layers, the mechanism is 
lost. The origin of the charge fluctuation effect will now be discussed in detail. 

We consider the Wannier representation of orbitals localized in a layer. The orbital wave- 
function has a representation 4>ji{r) without an explicit spin dependence. Actually, we can 
use a proper unitary transformation, e.g., the Wannier transformation, to create (pji{v) from 
Xm,'k,kz{^)- Here, I denotes an index of a layer and j represents a set of indexes (r?^, k). By 
associating creation and annihilation operators, and Cji are defined by the Canonical anti- 
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commutation relation, 

(^j'l'a'} = Sjj'6i^irda,a'- (21) 

As a relevant perturbation for the GGA band structure, we consider scattering processes 
coming from the charge fluctuation in (^iVxJ*). In the expression of the Y and Z operators, 
we have a pair of field operators, namely, V'o-(r) and V'o-Cr). These operators are expanded in 
the localized orbitals (l)ji{r) and ct^^, cjia- Thus, we notice that we have a double summation 
in the definitions of 1^^"^ and Z^'"\'^^) One is for the conserved quantities and the other is 
with respect to the orbitals in the Wannier form. Here, the center for the Coulomb expansion 
formula is not necessarily identical to a Wannier center. Depending on the symmetry of an 
electron pair in both the initial and final states, the proper selection of n = (p, L, M) is given 
to minimize the energy using DFVT. We can then identify relevant fluctuation terms, in which 

" in) 

Y-^ ' may connect a localized orbital in the low-energy bands at the Fermi energy with other 
semilocalized orbitals, leading to the pair hopping between Wannier centers. 

When a superconducting order parameter A;(r, r') = 4>il^)'\^i'l^'){c.ji'^Cjii\) is ex- 

pected to be finite, each term in can be re-expressed as below. As an example, we 

consider the expression 



Iihl3l4,,jij2j3j4 cr,a' 
~ E ^^9n;llll,jij2j3j4{- '^lila^jAlTCj^iaiCj^ia' ■) 

h=l2=l3=l4=l,jlj2j3j4, (T,a' 

+ E E^";^"'^'JlJ2-?3i4 {(4i;a42i<T')('^i3i'<T'Ci4i'a) 

h=l2=l¥'h=lA=l' ,31323334 <T,(J' 

+ E '^^dn^hhhhdihhhi'- ^jiha^34l4(T^j2l2a'^33h(r' ') 

h^h or 13^^14,31323334 o-,o-' 

- x\-\[y;:-^)'){y;:-^). (22) 

Here, the coeflficient 9n;hl2l3l4,ji323334 is given by 

9n;hl2l3i4,h323334 = ^ j J d^r'^*^i^{r) (#pLM(r) + $pLM(r)) ^j^ir) 

X cP*,i,{r'){%LM{r')+%LM{r'))cP,,i,{r'). (23) 
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A similar expansion is also given for terms with the operators. Thus, the effective two- 



body Haniiltonian Ti^ is divided into three parts: (1) the intralayer two-body Haniiltonian 
^f.intra (^ith gn;iiiLjij2r.m) ^ (2) the interlayer pair hopping Hamiltonian i^^ter-pair (with 
9n;m'i',h32hH)^ and (3) other 2-body terms Hf j^^ter-res- 

The first category "Hf^^^^^ contains terms interpreted as the on-site Hubbard repulsion. 
However, to obtain an explicit form, we need to introduce another Wannier transformation to 
have Wannier orbitals localized around a Wannier center in a layer. This process might be 
difficult, since we need to introduce unitary transformation in a rather wide energy window 
over a few eV. In our discussion, this process is not required to derive an effective low-energy 
model of high-temperature superconductors. In H^jjj^j.j^, we also have intralayer exchange in- 
teraction, intralayer off-site repulsion, and intralayer pair hopping. These effective interactions 
may induce the spin fluctuation effect as well as the charge fluctuation effect. For a general 
discussion, we do not specify the detailed form oi'H'j-^^^^^^ thereby allowing a BCS model Hamil- 
tonian for the electron-phonon mechanism, a correlated two-dimensional electron model, and 
a model considering both effects. In our pair-hopping mechanism, 'hL^ inter- -nair pl^ys a relevant 
role. 

(n) 

Note that is finite for a selected set of p, L, and M values. The second category 

inter-pair derived from a representation as a mean- field term plus the fluctuation term us- 
ing the superconducting order as in eq. (22). The third category 'H^ij^^.^j.^^.^g contains interlayer 
exchange interaction, inter-layer correlated hopping terms and interlayer diagonal charge fluc- 
tuation. We express the interlayer two-body Hamiltonian as = 'Hf^^^^^^_p^^^+Hf^^^^^j._j.^. 

The effective model may be written in a second quantized form as 

^eflf ^intra ^i.intra ^inter ~^ ^i,inter ~^ ^i,counter > (^^) 



'''intra 



'''inter 



EEE*5?'H4^^""^ + H.c.}, (26) 

where "T^intra the intralayer single-body Hamiltonian and T^i^ter the interlayer single-body 
Hamiltonian. The index / specifies an Z-th layer and j represents a j-th orbital in the l-th. layer. 
Note again that the j-th orbital is a Wannier representation made by one-dimensional Prourier 
transformation from the Bloch waves. Thus, we can always say that the interlayer single-body 
process for the Wannier states around the Fermi level is negligible for iron-arsenide supercon- 
ductors. At energy levels above the Fermi level, however, we also have finite bandwidths for 
the GGA band structure in the c direction. This picture is very important for the discussion 
below. 

We now construct a standard model of high-temperature superconductivity. Two steps are 
necessary. First is the construction of a mean-field description and second is the derivation 
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of the many-body effective Hamiltonian. We consider the effective intralayer Hamiltonians 
"^intra intra' which induce supcrconducting fluctuation. However, owing to their explicit 

two-dimensionality, Ti-l^^^^^.^ + Ti-^^^^j.^ cannot induce bulk superconductivity by itself. We then 
introduce the mean-field description of a two-dimensional superconducting state on a layer in a 
self-consistent field of other layers. Here, let us consider singlet superconductivity. The order 
parameter Ai{r,r') = 't'jl{'^)4'j'l{'^'){cmCjii\^) can have a finite value around the l-th. 

layer. One important point is that Ai(r, r') may change its phase as ( — 1)', but the following 
theory also allows a constant phase factor for all layers. The final superconducting phase 
should be determined by a variational determination method. 

We may consider two-body fluctuation, which induces pair-hopping processes between 
neighboring layers. The process comes from T~(-'j\^iQr-psXr- This term results in a pair field 
Hamiltonian for the layer. 




+^pLM(r)* X $pLM(r')) ^^iKr)0i'/(rOciiTCj'a + H.C.} 

= 5Z(^jV'^J'TC/U + H.c.}. (27) 
ji' 

(n) 

We have introduced an effective coupling constant ' for each scattering channel described 

by : {y^^^^ Y^^^ — [z\^^^ z\^^ :. It can be derived from the Coulomb kernel — how- 
ever, in the multireference density functional theory, effective coupling can be optimized. By 
applying the fluctuation reference method, we should determine the parameter for reproduc- 
ing another precise calculation. Or the effective Hamiltonian can be determined using DFVT, 
which will be addressed in the next section. A relevant point in our discussion is that the 
model is derived using the multireference density functional theory. 

The above derivation of eq. (27) is given by direct pair hopping. This Coulomb off-diagonal 
element is negligible for LaFeAsOi_j;Fj;, since the neighboring two iron layers are widely 
separated by a La2(Oi_a;F-r)2 layer. However, we have a pair-tunneling process across the 
insulating layer. We call it the super pair tunneling. To be precise, we show the construction 
step of the effective Hamiltonian for doped LaFeAsO. We can perform the non-magnetic GGA 
calculation of, e.g., LaFeAsO0.875F0.125 using a super cell with an optimized atomic position. 
The GGA band structure reveals the appearance of well-locaUzed 3d bands of iron at the 
Fermi level. Both electron and hole pockets are created from the localized 3d orbitals. These 
center bands have a clear two-dimensionality. Above these bands, at approximately 3 ~ 4eV 
higher than the Fermi level, we have delocalized bands that consist of non-s orbitals at La sites 
and void sites between a Fe2As2 layer and a La202 layer. These extended bands are formed by 
the hybridization between these high energy levels and iron 3d orbitals, so that localized 3d 
orbitals connect to higher levels using finite matrix elements by two-body Coulomb scattering 
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processes. With the help of these orbitals, an indirect pair-hopping process from a Fe2As2 layer 
to the next layer by Coulomb off-diagonal elements is allowed. This second-order perturbation 
process is relevant. Its effective form finally becomes the same as inter -pair' replace the 
interaction kernel with the effective one. We could also have a higher-order contribution from 
other terms in Ti-f for the pair-field Hamiltonian. However, it is important to note that 

inter always necessary to have the energy reduction in the superconducting state, since 
"^inter negligible at the Fermi energy and since Ti-f^^^^ is diagonal in the band index. 

For 11 compounds, we may consider only orbitals in a Fe2Se2 layer or a Fe2Te2 layer. In 
these systems, a direct pair hopping from one layer to the next layer is possible via Coulomb 
repulsion. Thus, we have two different categories of the pair-hopping mechanism. The first 
is the direct pair hopping and the second is the indirect super pair tunneling. Both of the 
processes require a finite amplitude for the pair hopping from a localized orbital to another 
well-defined orbital. 

The inclusion of the pair field necessarily results in a Josephson coupled superconducting 
state as a variational ground state. Its local wavefunction is given by the following effective 
Hamiltonian for the /-th layer: 

^Iff = E E ^1' Ha^j'i- + H.c.} + Hl^,,,^i + H^,, (28) 

(ii'> 

where the intralayer two-body Hamiltonian 'Hf^^^^^^i for the /-th layer can be either clectron- 
elcctron-intcraction originated, elcctron-phonon-interaction originated, or their combination. 
The pair-hopping processes producing can give energy gain to the Coulombic electron 
system, although a finite energy loss occurs when single-particle tunneling processes between 
layers induced by Hj^^gj, are terminated to have a variational state, even if they are negligible. 

In the mean-field description, we are able to obtain a mean-field solution, which have two 
order parameters: n(r) and Ai(r, r'). In the MR-DFT formalism, appearance of A;(r, r') affects 
the single-particle momentum distribution and its Fourier transform, i.e. n(r). Therefore, we 
can interpret that n(r) becomes A^-dependent. For the derivation, however, we need techniques 
for discretizing several continuous variables to have a tractable model in MR-DFT simulation. 
Here, we would rather move onto another effective theory to consider the superconducting 
phase derived from intcrlaycr pair-hopping processes. 

When a mean-field wavefunction of the layered material is obtained in the form 

i*(°))=nn«(-2.+-s.^L0io)> (29) 

I mk 

we can consider an explicit pair hopping. Here, ftj^j^^ = cj^kit'^m-kzj, '^ith cj^j^^^ = c^^^, where m 
is the band index and k is the two-dimensional wave vector. The real factor u'^-^^ and another 
complex factor, v^^^, satisfy (^^^-^ + l^nkl^ = 1- However, note that, even if we have a 
correlated superconducting state |^(°)) with an expression other than eq. (29), we can always 
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construct a Bogoliubov-Valatin transformation using the superconducting order parameter 
AKr,r'). 

Now we go to the final step in order to take the fluctuation efi'ect into account and 
to obtain the hidden order parameter and the energy gap. When a pair described by 
hops via the annihilation operation A;(r, r') = ^^jj') 4'jli^)4'j'li^')cji]Cjii^, the layer loses two 
electrons and the next layer obtains these electrons (see Fig. 1 (a)). Each layer should keep 
its charge neutrality, except for local charge fluctuation. Thus, if charge fluctuation effects 
are introduced into 1^''-''^), an effective action appears for this motion of pairs in the array 
of layers. (Fig. 1(b).) Then, we obtain a perturbed state l^'*-"'^-'), which is determined by the 
effective action of the pairs. We derive the effective action using only a unitary transformation 
without referring to these supposed state vectors. 

The charging effect should be of the same order of magnitude as the pair-hopping process. 
If the intralayer coherence is well kept, but if the interlayer single-electron hopping processes 
are not relevant, the screening effect expected for the system mainly occurs in the layer. The 
condition for Hj^j^j, is consistent with this picture. Then, the number of pairs allowed to hop 
at one time in a fluctuation process is restricted to be very small. 

To justify this discussion, we introduce the Josephson-Bardeen modification"^^' ^^-^ of the 
Bogoliubov-Valatin transformation: 

-t _ (Of _ (0 Q* 

-t _ (Oct _ (0 _ (? -t 

't _ (Oct j_ f (0 V _ <? -t 

where 5"; annihilates a coherent pair in the condensate of the Z-th layer and Sf creates one. It is 
not necessary to have the form of bl^^^ . Inserting the inverse transformation to a pair- hopping 
process, we find that 

'^Lk-t-pZT'^m-k-pU'^m'ki'TCm'-kZ'i 

_ / (0 , / (0 Y Q*- \( (0 "t _ (0 Q*- ^ 

~ l^"mk+p7emk+piT l^'"mk+p J '^l Temk+pU J y^'mk+p^emk+pll '"mk+p'^l 7emk+pZT J 

^ (('"mk+p)^7emk+piT'^Imk+pU + "mk+p (^mk+p) ~ '"•emk+pu) 



(0 (0 c*~ W ~ /o*\2- 

-""mk+p^mk+p-^i "emk+pZT " Vk+p K^l ) lemklllemkll 

(, (v) ' , in in c fi - \ c) f (nv o - 



,(0 
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Thus, we have a tunneUng process from a Cooper pair in a condensate in a layer to a pair of 
quasi-electrons in the next layer by the contribution (^^Sk+p)^"m'k^''k7Lk+pn7Lk+pU'^i' + 
H.c, which is found in the above expression. 

In a correlated quasi-electron system of the Z-th superconducting layer, the charge fluctu- 
ation mode in 'H^j^^,.g^ produces short-range repulsive terms for the quasi-electron state given 

by 

^emkZf ^^'^ ^emkll' ^^^^^ pair hopping occurs at a local position, the quasi-electrons in a 
pair inevitably raise their energy. To determine the localized nature of the hopping pair, we 
can do a rough estimation of the relative distance between two quasi-electrons in real space. 

d') w) 

For simplicity, we omit the radial dependences of ^i^/J^. and v!^/^, keeping only the energy 
dependences. Since we perform integrations in the k space, by considering a superconducting 
gap much smaller than the Kohn-Sham-band width of a few eV, we have the fohowing sim- 
plified expression for the creation operator of a hopping pair in an m-th band at an l-th. layer 
from the next Z'-th layer: 

J j cfkcfp (^^^c+p)^^^|7i'k'"i,'i7eTO,k+pit7lmk+pa 

~ // "^^^"^^^ ^i^rn,k+p - EF)diem',k - -^F)ik+piTik+pU 

- // d'kd^P ^(£m,k+p - EF)5{em',i. - ^F)4k+pn4-k-pu (31) 

Here, Ep is the Fermi energy of the Kohn-Sham band structure and we have omitted the kz 
dependence of Kohn-Sham energy owing to its two-dimensional nature. If we further consider 
a semimetallic band with the m-th conduction band, the above expression indicates that a 
contribution of doubly occupied Wannier states appears. In iron arsenide, this Wannier state 
should be in a localized 3d state at an iron site. 

The above-mentioned characteristic feature of the Hubbard-type correlated system is very 
important in considering layered superconductors with 3d local orbitals. The pair-hopping 
process from the condensate to a correlated quasi-particle state leads to the conclusion of a 
discretized energy spectrum as a function of the number of hopping processes at a time. In 
other words, depending on the number of quasi-particle pairs, npair, and mpair of pair holes 
in a layer, we obtain the energy contribution of (npair, ?TZpair) = C^'^pair + Uhm^aXv with the 
effective parameters ?7e > and ?7fe > 0. A simple correspondence of rzpair in the physics of 
the Hubbard model is shown by the number of doubly occupied sites. 

Let's therefore consider the situation given in Fig. 1(b), where only one pair is left from 
one layer to another layer. The bosonic nature of the pair allows us to write down an effective 
Hamiltonian in a Heisenberg spin system. Since we have a hopping pair or a vacancy at a 
layer, we have at least three states for each layer. We can assign these states to Sz = ±1 and 



^m'k 
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5*2 = of an artificial spin state and introduce an effective spin Hamiltonian. This minimal 
case is an 5 = 1 system, in which the ^-th layer may have one of these three states: Sz = 
and ±1. 

Now, the pair-hopping process is described by the xy-term in the Heisenbserg model. 
Neighboring charges with different signs will lower the energy, but two neighboring pairs with 
the same sign will raise the energy. This contribution is described by the anti-ferromagnetic 
z-term in the Heisenbserg exchange interaction. The charge neutrality discussed above also 
leads on-site anisotropy, i.e. the D term, to stabilize the Sz = state, which corresponds to 
a neutral layer. These contributions are described by an XXZ model with the D term. If a 
semimetallic band structure is the starting limit, and if the charge imbalance is minimized, the 
effective local charge neutrality is expressed by local energy enhancement, which is symmetric 
with quasi-particle pairs and pair holes, then a simple D term should appear. In general, the 
charge neutrality condition can be effectively expressed by the introduction of a D term. Thus, 
the minimal model is a one-dimensional anisotropic 5* = 1 Heisenberg antiferromagnetic spin 
chain. 

Wafhm = [JiSfSf+1 + SfSl,) + JzSfSU,] + D Y^iSff. (32) 
I I 

We know that there are three gapped phases in this model: the Neel phase, large-D phase, 
and Haldane phase.^^'^^^ The Neel phase corresponds to a pair- vacancy array in the present 
model. The state thus corresponds to a charge density wave state along the stacking direction. 
Thus, it is not relevant for the present consideration for the superconductor. The large-D phase 
corresponds to a Mott insulating phase, where the charge fluctuation effect is suppressed. In 
the limiting case in the large-D phase, a decoupled array of two-dimensional electron gas is 
realized owing to the suppression of the charge fluctuation. For a relevant contribution to sta- 
bilize a bulk superconductor, the Haldane phase is the necessary phase. In the Haldane phase, 
a valence-bond-solid state^^^ is realized with a broken hidden string order parameter.^^'^^) 
This gapped state with a hidden extra order contributes to the stabilization of bulk super- 
conductivity. The phase diagram of the antiferromagnetic S = 1 XXZ chain with uniaxial 
single-ion anisotropy is extensively studied.^^^^ 

The Haldane phase possesses an excitation gap. The lowest excitation with a total effective 
spin ^jQ^ = 1 corresponds to the creation of an extra Cooper pair in the bulk. In this state, the 
number of bosons actually increases by one. When the effective model has J = Jj = 1 with 
D = 0, the gap becomes A ~ 0.410479 J. There are continuous series of trials for determining 
Haldane gap. Recently, Ueda et al. have provided an estimation of lower and upper bounds 
and concluded that the gap is in [0.41047905,0.41047931].^) This study is performed by the 
combined use of the hyperbolic-deformation technique and sequence interval squeeze method. 
When interaction parameters are varied, the gap changes continuously in the Haldane phase. 
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Determining the value of the gap may be a simple test for checking the consistency of the 
present theory. If the gap is typically comparable to the transition temperature, we should 
roughly observe the pair hopping with J ~ /cbTc/0.410. If is approximately 100 K, J should 
be approximately 24meV. This value is within reasonable range for the two-body effective 
interaction, which is derived as a charge fluctuation. On the other hand, and D for diagonal 
elements of the effective model can have larger values without the electrostatic breakdown of 
vacuum or the insulating barrier layer between superconducting layers, e.g., Fe2As2. 

Here, note that an 5*2 = 1 state corresponds to any state with an extra pair in a layer, 
which may have multiple-colored states distinguished from each other. We only need to count 
the number of pairs that conies in or leaves a layer. If we consider states with many pairs 
coming in (or many vacancies leaving) a layer, we may utilize a higher-spin Heisenberg chain 
model. However, we know that Haldane gap decreases with increasing integer S. Thus, a high- 
temperature superconductor should be searched in a layered material, which can be mapped 
to an S" = 1 model. One might find that the above discussion is easily applied to triplet 
superconductors in layered materials, when inter-triplet-pair Coulomb fluctuation is taken 
into account. 

4. Density Functional Variational Theory 



Recently, the author has proposed a theory of the model space in the multireference 
density functional theory.^^'^^'^^) In this formulation, we use a new variational principle for 
the electron models defined by the density functional theory. A version of the density functional 
variational theory is given by the nequality 

< min jmin [*] + AEx,,e.,g. [*] | , (33) 

where Eo is the ground-state energy of the electron system, and * inserted in AExi,si,gi['^] is 
the minimizing of the functional Gxi,ei,gi[^]- The energy functional determining the model 
is given as 

+E,. [n*] + Eg. [^r] + J dVx;e^t(r)n*(r). (34) 

Here, E^^ [n] is an LDA energy functional, which may be a GGA energy functional, and Eg^ [^f] 
is a nonlocal correction parameter used in a standard DFT model. may be written 

by projection operators using a separable pseudo potential technique, an ultra-soft pseudo- 
potential technique, and a projector augmented wave technique. We introduced an energy 
difference functional as 

AEx„e„gA^] = ^ / d'^dh'j^^i'^l ■■ (^(r) - "*(r))(n(r') - n^(r')) : 1^-) 
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- E,^[n^]-Eg^[^]-mVx^^). (35) 
The proof of eq. (33) is easily realized by noting the next inequality 

- mil {(*'|t+l>5|*')+B„|9/']} 

= YjdxJ ciVdV^^«| : (n(r) - n{r)){n{r') - n(r')) : |0 
= min (*'|f + t>ee|^') - - / d^rd^r'^^^P^ 

«"^n* 2 J k - I" I 

2 J |r — r'l 

- E,,[n^]-Eg,[^']-{^\f + V^''^^\^) 
/TiT> I TV /" j3 j3 /"*(r)'^*(r') 

= (^l^eel^') - Y y d^^ ^ 

= ^Ex„e„g,m. (36) 

Here, we used 

min {(*'|f + V;fj|*') = (*|f + V;f^|*) (37) 

since the minimizing |^') of the above expression for the charge density n^(r) is obtained by 
minimizing a functional, 

+^,Jn>i,] + E^J*'] + y d^rvea^tir)n^ir). (38) 
which is given by |^'). We again find importance of the self-consistency in the minimizing 

process of Gxi,ei,gi[^]- 

Thus, 

£^0 < min \imnGxueu9i[^] + AEx,,ei,gA^]] 

< min \mmGx„e„9iM + AEx„e„g,m \ . 

Thus, we can start from a known LDA functional to construct a variational model of the 
electron system. The fluctuation term Vx^ is formed by static two-body correlation functions. 
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which yield two-body effective interactions in the MR-DFT model. Thus, we have a firm 
ground for a beyond-LDA approach, considering that relevant fluctuation modes are inserted 
in the variational model. We may utilize GGA in place of LDA under the assumption that 
the exchange-correlation potential is always given in the simulation process. 

In the density functional variational theory (DFVT), a simulation for determining the 
fluctuation term (^'|Vx,. |^) is given, when a differentiable ii'ejn] and/or an explicit Eg-l"^] are 
prepared. Then, we always have a defined one-body part of the effective Hamiltonian in a self- 
consistent determination process of the self-consistent solution of the model. Thus, an LDA 
or GGA solution can be used to construct a correlated electron model of superconductivity. 
The model has effective two-body interaction terms. According to the model Hamiltonian, 
we can apply any appropriate solver for an effective many-body problem. Self-consistency is 
imposed by calculating the charge density, which redefines the exchange-correlation potential. 
Two-body processes and their interaction parameters are determined, so that they reduce the 
variational energy of the Coulomb system. The variational energy is given by evaluating all the 
terms in eqs. (33) and (35). At present, realistic determination processes are computationally 
demanding, but a preliminary simulation of Sr2Cu03'^^^ indicates that an optimization pro- 
cess indeed works for determining an effective interaction parameter. In this one-dimensional 
copper oxide with the (P configuration at each Cu, the on-site Hubbard interaction J7 > is 
determined by searching the minimum variational energy. 

5. Conclusions 

We proposed a pair-hopping mechanism expected in a layered superconductor on the 
basis of MR-DFT. The definition of the effective model for a superconductor is given from 
the first-principles method. Determination techniques of the model is given by DFVT. The 
derived effective pair-hopping model suggests the realization of a valence-bond-solid state of 
an effective S = 1 Heisenberg anti-ferromagnetic spin chain. A bulk superconducting phase 
becomes gapped in the charge fluctuation mode. 

Sufficient conditions for the pair-hopping mechanism are summarized as follows: (i) There 
is a two-body pair hopping process between layers, which may be direct or indirect, (ii) Inter- 
layer single-particle tunneling is negligible compared with the two-body pair hopping between 
layers, (iii) A correlation effect in a layer suppresses multiple-pair hopping at the same time, 
keeping the local charge neutrality and allowing minimum charge fluctuation in a uniform 
bulk superconducting state. These conditions allow us to have a high-temperature supercon- 
ductor. The realization of the mechanism in iron arsenides is expected because we obtain 
(1) experimental observation of high-temperature superconductivity, (2) the two-dimensional 
electronic state given by GGA, and (3) the present formulation of the pair-hopping mechanism 
in MR-DFT. DFVT suggests that, if competing diagonal orders in magnetic and non- magnetic 
channels are not comparable in variational energy, the superconducting state is selected. For 
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the ofF-diagonal superconducting order, there is energy reduction in Coulomb energy, because 
the pair-hopping mechanism selects the unique ground state with the Haldane gap. 

The importance of the interlayer pair tunneling process has often been stressed for cuprate 
high-temperature superconductors. "^^"^^^ In our DFVT, we can also derive a strategy for en- 
hancing the stability of the superconductivity using our microscopic effective model. The 
Coulomb-originated pair hopping is derived via the charge fluctuation modes. One possible 
form is given by eq. (5). Here, we need to consider the sign of the superconducting order 
parameter. We have two contributions, namely, a ^-process and a Z-process, due to the ap- 
pearance of two operators. The important point is that the sign of these processes are different. 
Depending on the two-dimensional superconducting order parameter, one of them can be ef- 
fective for direct pair hopping. As for the super pair tunneling mechanism, the combination of 
the y-process and Z-process may appear in the higher-order pair hopping process from one 
layer to another layer. 

In Fig. 2(a), we show a possible scattering process from one layer to another layer. By 
using a l^-process, two branches of the gap function, i.e., an m-th band with a positive sign 
and another m'-th band with a negative sign may be stabilized. In a Z-process, the final pair 
potential does not need to change the sign from the initial one. The momentum conserva- 
tion at a scattering center holds. According to this rule, we can find out relevant scattering 
processes depending on the geometry of the Fermi surfaces and the sign of the local order 
parameter ^jj'i- Two examples shown in Fig. 2 are for (b) an extended s-wavc state modeling 
iron arsenides and (c) a d-wave state for cupratcs. The strength of these processes is depen- 
dent on the orbitals and effective coupling strength, and thus on the material striicture 
considered. Furthermore, we need a strong superconducting fluctuation in a two-dimensional 
layer. For the enhancement of the intralayer fluctuation, we may rely on our knowledge of 
the spin-fluctuation mechanism of high-temperature superconductivity.^'^' Another known 
fact in the literature supporting the present pair-hopping mechanism is the electronic band 
structure calculations for optimally doped LaFeAsO,^'^-' which has a higher Tc than NaFeAs. 
The two-dimensionality of Fermi surfaces is confirmed, when the experimentally observed lat- 
tice structure is assumed in the simulation or when hole doping is assumed in a theoretically 
determined lattice structure. A perfect two-dimensionality of the LDA or GGA band struc- 
tures means a strong suppression of interlayer single-particle hopping processes as well as the 
appearance of correlation effects in each layer. The indirect pair tunneling mechanism further 
suggests the co-existence of magnetic order in an independent part of the layered supercon- 
ducting system. If the magnetic structure only acts as a medium supporting the pair hopping 
processes, the gap formation in a stack of local superconducting two-dimensional electron 
systems is expected. 

A simple comment on another attempt^^^ to elucidate a superconductivity is given for 
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pedagogical reason. The present pair hopping mechanism is aUowcd only for the interlayer pair 
scattering due to the Coulomb repulsion between electrons. If there is an effective attractive 
interaction at an attractive center between layers, the transition temperature vanishes, since 
the pair hopping is blocked. The existence of a on-site static attractive interaction coming 
from the Coulomb repulsion has been already disproved. 

Finally, three possible comments on real superconductors are given as follows. The origin 
of inter-layer scattering processes relevant to bulk superconductivity may be the electron- 
phonon interaction as well as the intra-layer effective attraction. In the case of MgB2,^'^-' the 
pair-hopping mechanism can give some amount of stabilization through the appearance of a 
hidden order parameter. Since we need to specify details of scattering processes that stabilize 
the Coulonibic electron system, further theoretical investigation might be necessary. 

For the realization of the effective S" = 1 Heisenberg model, we conjectured that a 
semimetallic band structure is favorable. Pseudo-electron-hole symmetry corresponds to the 
XXZ Heisenberg chain with the D term. To have the most stable Haldane gap, we also need 
to perform numerical simulation of a generalized Heisenberg chain model. Although a detailed 
discussion on the stability of the multiple-ordered state proposed on the basis of the present 
pair-hopping mechanism is required to determine the best condition by first-principles simu- 
lation, we expect to obtain a reliable estimator of Tc in the near future, since the numerical 
accuracy for the determination of the Haldane gap is now going beyond single precision. '^^^ 

In several cupratc superconductors, the realization of microscopic Joscphson junction ar- 
rays is known to be related to the Joscphson plasma. "''^"^^^ For the microscopic analysis of 
this phenomenon, the determination of the sign and structure of the superconducting order 
parameter is important. A simulation based on DFVT is expected to solve this problem, too. 
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Fig. 1. Minimal model of the high-temperature superconductivity in a layered material, (a) In the 
layered superconductor, a Cooper pair may hop via a two-body pair-hopping process into the 
next layer, but single-electron tunneling processes are suppressed, (b) Extra pairs and resulting 
vacancies may be created by the hopping process with charge fluctuation, (c) Neighboring pairs 
or neighboring vacancy layers are energetically forbidden. Thus, we have a hidden string order in 
the Haldane phase. 
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Fig. 2. (a) Effective interaction process for the pair hopping from one layer to another layer. In the 
first /-th layer, a Cooper pair experiencing a positive pair potential is scattered into the next 
Z'-th layer by a F-process, where the final state may be affected by a negative pair potential. In 
the scattering process, an electron with (m,k, Z,(7) of a pair is scattered into (m', k', /', cr) with 
a momentum shift and another electron (m, — k, I, —a) of the pair obtains an inverse momentum 
shift to be (m', — k', V , ~a). (b) Schematic viewgraph of the superconducting order parameter in 
the two-dimensional first Brillouin-zone of a model two-band superconducting state. If the Fermi 
surfaces are separated into a center one and another one around the M point, an extended s-wave 
state may be created. The interlayer pair hopping from the /-th center Fermi surface to the I'-th 
Fermi surface around M may occur, (c) A schematic viewgraph representing the sign of the order 
parameter for a model superconducting state with a d-wave state. 
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